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On the Matrix which Represents a Vector, 

By C. H. Chapman. 



The fundamental idea of this work is emphasized by Buchheim as follows : 
"It is, or ought to be, well known that the linear and vector function of a 
vector is simply the matrix of the third order,"* meaning, of course, that the 
matrix is the operator which transforms p into $p (Tait's Quaternions, p. 98. I 
shall refer to Tait's work from time to time simply as "Tait," the edition being 
that of 1867). It is merely a very special application of this remark to observe 
that any vector whatever whose rectangular components have the lengths x, y, z 
may be derived from that of which the components are of length 1 , 1 , 1 by the 
operation of the matrix 

( x ) 

y 
2 

and it follows that the study of vectors may be made to depend upon that of 
these matrices. The object of the present paper is to carry out and illustrate 
this study in some detail. 

In the course of the work the matrix is freely spoken of as the vector by an 
allowable looseness of language, since no confusion can result; but it is not 
assumed that they are the same — in fact I wish to guard most carefully against 
that assumption. 

The algebra of these matrices is the same as that of scalar quantities, inas- 
much as, owing to their very special form, they are commutative with one 
another in multiplication ; for this reason I have found it necessary to employ 
only a few very elementary properties of matrices. In fact Oayley's celebrated 
memoir (Phil. Trans. 1858) contains nearly all the theorems which are here 

* Proo. London Math. Soc, Vol. XVI, p. 63 ; Clifford, Dynamic, p. 186. 
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made use of, although I have referred by preference to the papers by Sylvester 
(Am. Journ. Math., Vol. VI) and Taber (Am. Jburn. Math., Vol. XII) as being 
more recently in the hands of all readers. The learned article of Taber contains 
a statement of pretty nearly the entire theory of matrices, at least from his 
point of view. The only new symbol here introduced is that of a circular sub- 
stitution, cy or cy' , which has proved very interesting and useful ; for the rest, 
the symbols S and V, not as selective but merely appellative, but with their 
geometrical significance essentially unchanged, have been found sufficient. 

Many, if not all, the advantages which result from the use of quaternions 
will be found attained or attainable in this paper without their aid by the use of 
these matrices ; the equations and formulas are sometimes identical in appear- 
ance with those of the quaternion calculus, and the resemblance has been syste- 
matically cherished — partly to aid in making comparisons, partly on account of 
the extreme elegance attained by the quaternion calculus in the hands of Hamil- 
ton and Tait. 

The unpublished Vector Analysis of J. Willard Glbbs is briefly referred to 
in this article as V. A. Gt. 

1. — Properties of the Matrix. 



p is referred to its axes,* which will be assumed 

as a rectangular system and, for distinction, designated as the axes of i, j and 7s 
respectively. The matrix p denotes that substitution which transforms the 
vector whose components along the three axes have the lengths 1,1,1 into that 
whose components are of lengths x, y, z respectively. The matrix therefore 
defines without ambiguity the length and direction of this vector, and I shall 
say that the matrix is the analytical expression or representation of the vector. 



( 


X 





o) 


The matrix 





y 













z 



In particular the matrix 



1 











1 











1 



represents a vector whose length is \/3 



and which makes equal angles with the positive directions of the axes. In the 



* Taber, Am. Journ. Math., Vol. XII, p. 360. 
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vector theory this matrix plays the part of unity in the scalar and will be spoken 
of as vector unity and denoted by v.* 

For brevity, matrices of the kind here considered will be denoted by Greek 
letters, and the same letter will uniformly denote both the matrix and the vector 
represented by it. If 

( a ) (a' 



a = 



then 



a + £ = 



b 

c 

( a -f- o! 


, 



and /3 = 

, 

b + V, 

, 



) 
V 









c + d 



)t 



Whence it is evident that the sum of the matrices represents the sum of the 
vectors. J 

The product of the matrices a and (3 is 

( aa! ) 



a(3 = 



bV 





GG 



which is evidently the same as fia, and, for the purposes of this paper, the 
product of the vectors a and /? will be so defined that it may be represented by 
the matrix a(3; that is, this product is a new vector whose component lengths 
are the products of those of the factors. The writer is well aware of his pre- 
sumption in thus departing from illustrious usage, § but inasmuch as he thus 
attains analytical consistency, simplicity of notation, and preservation of the 
associative law, both of which latter are lacking in the ingenious work of J. Wil- 
lard Gibbs, and at the same time retains the freedom of commutative multipli- 
cation, the lack of which makes Hamilton's system at once so significant, so 
fascinating, and so difficult for ordinary students to master and apply, he has 
ventured, in a tentative way, to make the departure. Unless all the constituents 
of the matrix except those of the principal diagonal are zeros, a/? is not in 
general the same as /?a, so that I by no means imply by this that linear vector 
functions are in general commutative in multiplication. 



* Sylvester, Am. Journ. Math., Vol. VI, pp. 274 and 275. 
t Sylvester, Am. Journ. Math., Vol. VI, pp. 274 and 275. 
% V. A. G., p. 4. § Tait, p. 44 S. ; V. A. G., p. 5. 
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Since division of matrices is defined by the equation 



we see at once that the quotient of the matrix a by /3 is the matrix 



( a 



^-0 

a! 





) 



b_ 

V 





c 



which will be taken to represent the quotient of the vector a by the vector /3 . 
From this it follows that the division is not ambiguous and that 

It is to be observed that this division is indeterminate or gives a quotient with 
an infinite constituent in cases quite analogous to those of scalar division. f 

Without using Sylvester's Theorem, J it is evident from the rules for multi- 
plication and division that any function /(p) which can be expanded in positive 
or negative powers of p will be denoted by the matrix 

(/(*) o o ) 

o f(y) o 

/(a) 



The matrices 



(1 















> ( 


0" 


> ( 








> 


10 


) 















1 



will be denoted by i, j, h respectively. For their multiplication we have 



i 2 — i t j z =y , h % =7e, ij = ih =jh = 









o, 





















= ijh. 



* Sylvester, Am. Journ. Math., Vol. VI; p. $76. 

t V. A. G., p. 47, Art. 128; Sylvester, Am. Journ. Math., Vol. VI, p. 274. 

% Quoted and verified by Taber, p. 378; Buchheim, Proc. Lond. Math. Soc, Vol. 16, p. 81. 
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Also ip = xi ; jp = yj; &p = zJe; 

whence (i + j + &) p = vp = p = ip +/p + hp ; 

being a formula for decomposing a vector into three rectangular components. 
The matrices 

( 1 0") (000) ( 1 ) 



10 





10 
1 




1 



= I,J,K, 



have the properties that 



P = I, J* = J, K* = K, 
IJ=JI=j- IK=i; JK=h, 
il=i; iJ= 0; iK=zi, 
jI=j;jJ=j;jK=iQ, 
M=0; JcJ=Jc; hK=Tc. 

The vectors i,j, Te; I, J, K cannot be used as divisors. 

2. — The Operators cy and cy' and the Symbols S and V. 
(x ) 



If p denote the matrix 



, or briefly (x,y,z), then shall cyp 



y 

z 

denote (y, z, x) and cy*p, (z, x, y); cy being an abbreviation for the word cyclic 
and denoting evidently a circular substitution among the constituents of p.* 
Clearly also c „3 p __ _ ^ 

Again, if a , (3 , y are any three rectangular unit vectors and p = xa-\-y(3-\- zy , 
cy' shall denote the substitution 

(<*> 0i V\ 7, «i /?)> 

so that cy'p = xy -{- ya -\- zfi ; cy'*p = x@ + yy + za, cy''p = p. 

Denoting a + /3 + y by »', cy and cy' are operators which rotate p in a negative 

direction about the axes v and xt respectively, through an angle equal to 



In 



while the rotation performed by cy % and cy'* is 



4n 



*Cf. Taber, Am. Journ. Math., Vol. XIII, p. 163. 
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If a = (a, b, c) and /? = (a', V, c') are any two matrices of the kind in ques- 
tion, then a/3 = (aa 1 , W, cd) and 

S.a(3 = ad + W + cc! 

will be called the scalar of a(3 . It is the product of the lengths of the vectors 
into the cosine of their included angle, and has thus its usual geometrical mean- 
ing. Its practical use is the same as in quaternions. Since in any product the 
vector v may be suppressed, the expression S.vp = x + y + z will often be 



written simply 


S.p. 


Obviously 
Also 




cySp = Scyp and cyS.a(3 = Scy (a(3). 
c#(«, P,y,8 ) = cyacyficyy . . . 


and 




a cya 
Cy P ~ cy(3 * 



While cy (acy(3) = cyacy^ft and S. acy % (3 = S. {3cya, 

from the definitions of the symbols. We recall also that if S.a(3 = 0, a and (3 

are perpendicular to each other, and note that the product of any two vectors 

perpendicular to each other is a third vector perpendicular to vector unity. 

In particular, 

cyi—Jc , cyk=j, cyj —i, 

cyl— K, cyK= J, cyj= I. 

Let us now form the vector 

8 = cyacy 2 (3 — cy(3cy 2 a. 

We have by the above formulas 

S.a (cyacy % (3 — cy^cy^a) — 8. acyacy % (3 — S. acyacy^fi = ; 

and likewise S. (3 (cya cy 2 (3 — cy(3 cy 2 a) = . 

Hence $ is perpendicular to a and {3 . Its direction is such that a positive rota- 
tion about it as an axis turns a to /3 , and its length is the product of the lengths 
of a and /3 into the sine of their included angle ; it is therefore the vector which 
Hamilton denotes by Va(3 ; but since we do not derive it from a and {3 by multi- 
plication, it will be denoted in this paper by V(a, /?), or — V{(3, a). 

We have V(i,j) = Jc; V(j, k) = i, V(k, i) =j. 
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Also, if a, /?, y are any rectangular system of unit vectors, 

V(a,(3) = y, F(y f a) = 0; V(fi, y) = a. 
We observe now that taking a + /? + y = v' for the axis of c?/, then 

cy'a = y, ciffi — a, cy'y = {3; 

and that if p = tea + «//? + ay , then 

S.tfp== x -{■ y -{■ z. 

Finally, however p may be expressed, the square of its length is 

It will now be clear that 

S.aVtf, y) = S.(3V(y, a)=S.yV(a, /?),* 

whatever vectors a, {3 and y may be; and that if S.aV(@, y) = 0, the three 
vectors lie in the same plane. 

The coefficients x, y and z in the expression 

8 = xa + yfi + zy, 

where a, /?, y do not lie in the same plane, can now be determined. Multiply- 
ing by V(a, /?) and taking the scalars, we find 

S.8V(a,P) 
Z -S.yV(a,@) {y) 

and similarly 

_ S.&V(P,y) __ S.8V(y, a) 

X ~ S.aV(J3,y)' y -S.(3V(y, a)' 
Whence 

&S.aV(P, y) = aS.8V(P, y) + pS.$V(y, a) + yS.8V(a, (3).f (1) 

This expression is scarcely less simple than the corresponding one in quaternions, 
from which it differs in form only by the notation V(a, fi) instead of V.a(3. 
In the quaternion formula S.(aV.@y) the V may be omitted without changing 
the meaning, but from /S.aV(ft, y) it cannot. 
It is also possible to expand 8 in the form 

S = xV(y, a) + yV(a, P) + zV(@, y). 

* Cf. Tait, pp. 56 and 65. t Tait, p. 57. 

47 



370 Chapman: On the Matrix which Represents a Vector. 

We have 

__ S.aS __ S.y8 S.(38 

Z ~ S.aV(P, r )'> y- S.yV{a, f) ; X ~ S.pVfr, a) ' 
Hence 

$S.aV(P,y)= V(y,a) 8.08+ V(a, (3) S.yS + V(P,y)S.a8. (2) 

The following formulas which follow directly from the definitions of the 
symbols cy and S are useful in what follows : 

vS.p = p + cyp + cy 2 p, \ 

vS.p* = p 3 + cyp* + ctfp\ I (3) 

whence v (S*. p — iS. p 2 ) = 2vS. pcyp . ) 

Also v?S. v'p = p + cy'p + cy'*p . (4) 

The cosine of the angle which p makes with v is 

VSiS.p"' 

and the cosine of the angle between p and v' is 

S.v'p 
V3£y* 

The components of p parallel and perpendicular to v are respectively 

^vS.p and p — ^v./Sp; (5) 

while the same with respect to v' are 

^v' /S.v'p and p — \v' S.v'p. (6) 

3. — The Symbols cy 2 * and cy®*. 

Since cy 3 and cy 1 ' perform respectively negative rotations equal to 2n about 
their respective axes, it is very natural to indicate corresponding rotations 

— — 2/7t 

through the angle 3 by cy 2 * and cy' 2 *. When S- = , 



3* 3& 

cy 2 * = cy , cy 12 * = cy'; 
when $• = — — , cy 2 * = cy 2 ; cy 12 * = cy'*. 

o 
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I propose to find a matrix (i such that 



3* 



Evidently for $• = 



2?t 



^ 



and for 3 = — , 
3 ' 



^ 



w- 


= cy^'p . 


(JL 

X 






















X 

z 


(JL 





) 





IT 











2 



(7) 



(8) 



We observe that p and ftp have the same length, and make the same angles with 
v; hence we obtain the equations 

S.pp =S.p; (9) 

S.tfp % = S.f. (10) 

To obtain the third necessary scalar equation, we note that §■ is the angle 
between ^ (3p — vS. p) and ^ (S(ip — vS.fip) , the vector projections of p and ftp 
respectively perpendicular to v . Hence 

rn ^_ iS.(3p-vS.p)(S W -vS.(ip) _ 3S. W i -S.pS.(ip 
cos '" %S.(5p — vS.pf " _ 3/S.p 3 -^.p " 

Whence 

ZS.pf = (3S.p 2 - S*.p) cosS + S 2 . p. (11) 

Instead of solving these three equations directly for /i , let us approach the result 
indirectly by aid of a vector v which rotates (ip through an angle ^>, leading to 

the equations S.vpp = iS.pp = S.p, (12) 

S. vYp* = S.fiY - S. p\ (13) 

SS.vfiy =(3S.p % — JS\p)cos<p + l S\p. (14) 

Again we notice that v(i is a vector which rotates p about v through the angle 
$ + <?> without altering its length, and thus deduce the additional equations 

ZS.vtf = (SS.p 2 - S*.p) cos (S + <?>) + S*p. (15) 
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Now in general the vectors (ip, (ip 2 , and (i 2 p 2 do not lie in the same plane, for we 

have 

vS. (j?p 2 V (ftp , fip 2 ) = vS. ^cypcy'^cypcy'p (cy 2 p — cyp) , (16) 

remembering that cyp 2 = (cyp) 2 and cy 2 p 2 = (cy'p) 2 . And if we take into account 
the fact that 

vS. pcypcy 2 pa = pcypcy 2 pS. a (17) 

we may write equation (16) 

vS. fa* V (ftp , (ip 2 ) = ncyiicy^pcypcytpS. pp (cfp — cyp) . (18) 

If S.pp (cy 2 p — cyp) = , then (i satisfies a fourth scalar equation besides (9), 
(10) and (11). This can happen, for a given value of p, only in very special 
cases, and one such case is {i = v . 

In general, then, by aid of equation (2), we may write 

v.S.[ipV({ip\ nY) 

= VQiy, up) s. vrf + v dip, ^p 2 ) s. v/xy + v (^ P 2 , fiy) s.v^. (19) 

In simplifying the formula (18) it will be useful to note that 

V (8a, 8(3) = cy8cf8 . V(a , /?) . (20) 

We may see this to be true since 

V(8a, hp) - cy (8a) cy 2 (5/3) — cy (5/3) cy 2 (8a) 

= cy8cy 2 8 (cyacy 2 fi — cyficy 2 a) = cyScifS .V(a, /3) . 

Using this formula, equation (19) becomes 

vS-upcy (fip 2 ) ctf (pp 2 ) V(v , n) = cy (fip) cy 2 (fip) V((ip, v) S.Vfip 2 

+ cy (jip) cy 2 (up) V(v , p) S. v{i 2 p 2 + cy (jip 2 ) cy 2 (ftp 2 ) V(v , p) S. vfip . (21) 

Using the formula (17) we see that this equation is divisible by the factor 
cy{ipcy*iip, giving for the quotient 

vfipS.cypcfpV^, (i) 

= V((ip, v) S.vpp 2 + V (v, p) S.v(i 2 p 2 + cypcy*pV(v, (i) S.vpp. (22) 

To complete the determination of v we shall assign to ^ a particular value 
such that knowing the effect of von /ip, we can readily find the vector which 

rotates p through the same angle. The vector (i = -^- satisfies equations (9), 

P 
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(10) and (11) but not (18), except in the special cases when p is the vector to 
the surface xy + yz + zx=ot? + tf + s 2 . 

We may therefore take u = -&* , and the corresponding value of 3- is -5- . With 

p - 6 

these values equation (22) becomes 

SvcypS.cypcfpV(v, 3£) = V{cyo, tf)[(3tf.p 2 — S*.p) cos (^- +<?>) + S\p~\ 

+ *>> P)[(3^.p 3 - £ 3 -p) cos^ + ^.p] + 3otfpo^pV(v, Wtys.p. (23) 

cypcy 2 p F(v , ^£) = cypcfp (%£ - ^£) = pcyp - cy 2 p 2 , 

and #. (pcyp — cy 2 p 2 ) = S.pcyp — S. p 2 ; 

since S.cyY = cifS.p* = #. p 3 . 

And from equations (3), 

£.pcyp=£(# 2 .p-,$.p 2 ). 

HenCe £.cypcy 2 pF(«, ^P)= i(£ 2 .p-3tf.p 2 ). 

Again, 

(T(cyp, ») + F(t>, p)] £ 2 .p + 3cypey 2 pF(V 32) tf.p 

= V(v, p - cyp) S\p + 3F(p, cyp) S.p 

= [V{p + cyp + cy 2 p, p — cyp) + 3F(p, cyp)] tf.p 

= IT(p, ^p) + Wp, p-cyp)] S.p=hv [S\p- 3£.p 2 ] S.p. 

Dividing out the common factor [/S^-p — 3S.p*] from equation (23) we have left 

ivcyp =: — T(cyp, «) cos (-^ + <M — "P"(^, p) cos$ -f- i v S.p, 
or finally, 

Svcyp =2V(v, cyp) cos (-— + <p) + 2V(p, v) cos$ + vS.p. (24) 

Changing cyp to p and consequently p to cy a p, and recalling that 

S. p = S. cyp = /S. cy 2 p , 
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we have, if we replace v by fi , 

-y + $) + 2F(c2/ 2 p, v) cos$ + vS.p. (25) 

Now, inasmuch as v was a vector which rotated cyp through the angle $ , it must 
be that p has the same effect on p . 

Equation (25) is homogeneous in the tensor of p and (i is therefore not a 
function of the length of p , but depends merely on its direction. 

As a partial verification of the correctness of this formula, I remark that 

when $ = it becomes 

3^p = 3p , or [i = v, 
as it should. 

Since (25) gives one and only one value of n for each value of $, no two 

being the same, as $ varies from to 2n, [ip will perform a complete rotation 

about the axis v and return to its original position. 

We may now write 

f* (26) 

cy*y = [if>, 

where ^ is a known vector. 

In a similar manner we may compute the vector factor which will produce 
a given rotation about any given axis. 

In fact, denoting any vector axis of length \ZSbya + /3 + y, where a, /S 
and y are three rectangular unit vectors so situated that the given axis is the 
diagonal of a cube of which a, /?, y are the edges, any vector whatever may be 
denoted by p = xa + y(3 + zy, (27) 

and the problem is to find a substitution p which rotates p through an angle $• 
about «'. 

The equations of condition are 

S.iiptt = S.fvT, (28) 

S.fiY = S.p z , ( 29 ) 
S.rf — i S*.vw = cos S (#.p 2 — i S'.vfp), 

or S. w* = | [cos a (3#. p 3 - S 2 . t/p) + S*. t/p] . (30) 

Equation (30) is obtained as follows : 

The component of p parallel to v' being tv 1 , we shall have S.tf (p — tv') = 0, 

whence tS.v 1 * = S.v'p, 

or t = ^S.v'p. 
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Hence the component of o perpendicular to v 1 is 

p — WS.t/p, (31) 

and that of [ip is [ip— ^v'S.tffip. (32) 

The vectors (31) and (32) have the same length. Observing that S.v** = 3, and 
that the scalar of their product divided by the square of the length of either 
gives the cosine of the angle between them, we have 

. _ S .iip 2 -^ S 2 .v'(ip 
G ° 8S >- s.p* — is*.i/p ' 

from which equation (30) follows at once. 

Proceeding as before we shall obtain three equations corresponding to (9), 
(10) and (11) respectively, 

S.vnp* = 1 [cos (3. + <p)(3S.p 2 — S\v'p) + # 2 .*/p] , (33) 

S. vfiy = i [cos $ (3S. p 2 — S^p) + S*. t/p] , (34) 

S.i/izvp=S.v>p. (35) 

Using the formula (2) to expand v in terms of ftp 2 , [i 2 p 2 , and i/pp, we obtain the 
equation 

vS.xttipV^p*, fa 2 ) 

= v((iy, fipv 1 ) s.v(ip* + VQipv 1 , [ip*) s.vfiy + v^p\ ^y) #.«vp, (36) 

which, by aid of (13), becomes 

vS.v'fipcy (^p 2 ) cy 2 Qip 2 ) V(v , (i) = cy ((ip) cy 2 (jip) V(pp , v') S.vpp 2 

+ cy dip) cy 2 (jip) V(vf, p) S. v(i 2 p 2 + cy (up*) cy* (jip 2 ) V(v, [i) S.v'vfip. (37) 

From this we may divide out the factor cy (pp) cy 2 (jip) , leaving 

v[ipS.v'cypcy*pV(v, (i) 

= V((ip, v') iS.vfip* + V(v', p) S.v[/p 2 + cypcy % p V (v , fi) S.v'vfip 
= $S».v>p.[V([ip,vi)+V(vi,p)] ' }-(38) 

+ £ (3#.p 2 — aS 2 .«'p) . [7 ([ip, v).cos (3 + <£) + V(v', p) cos $] 
+ S.v'p.cypcy 2 pV(v, (i) 

by aid of equations (33), (34) and (35). 
Equation (38) may be written 

vfipS. v'V(p, (ip) = £ S 2 . v'p . V (v/, p — ftp) 

+ i(3£.p 2 — S 2 .v'p)[V(np, *>') cos(3 + <?>) + V(v', p) cos<?>] + V(p, {ip)S.v'p. 
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If in this we assign to ft the value ?&-£ and to S- the corresponding value of — — - , 

9 3 

it becomes 

vcy'pS. v' V (p , cy'p) = J S 2 . v'p . V («/, p — cy'p) 

+ i (3#.p 2 — S z .v'p)\_V(cy'p, v') cos (^ + <?>) + V(J, p) cos^>] -(39) 

+ V(p, cy'p) S. v'p. 

The vector v turns cy'p through the angle ^> around the axis v' ; hence changing 
cy'p to p and therefore ci/'p to cy'p and p to cy'*p , we shall have the vector ft 
sought for : 

(ip/S.v , V(cy'*p, p) = £ S z .v*p. V(v', cy'*p — p) 

+ i (3tf.p a — S».v>p)[ V(p, v') cos (-— •+•<?>) + F« c y 2 p) cos <?>] - (40) 

+ V(cy'*p,p)S.<tfp. 

Now c«/' 2 p = x/3 + £/y -f- z<x and T 7 " (c«/' 2 p , p) 

= (xcyfi + ycyy + zcya)(xcy z a + ycy^fi -f- zcy 2 <y) 

— (xcy z /3 + ycy*y + zcy z a)(xcya + ycyfi + zcya) 
= {yz--x*)V{a, @) +(zx — y*)Vtf, y) + (xy'—z*)V(y, *) 
= (zx — y % )a + (xy — z 2 ) @ + (yz — x*)y, 

owing to the property of three rectangular unit vectors that 

V(a, 0) = y; F(0, y) = a; F(y , a) = /?.* 
We conclude that 

/S.v'V(cy' a p, p) = (zas + xy -{- «/s) — (as 2 -J- s/ 2 + a 2 ) = — i (3#.p 2 — jS 2 .v'p). 
Also 
F(«', c y 2 p — p) = V(v>, {z — x) a + (x — y) (3 + {y — z)y) 

= (cya + c«//3 -J- cyy)((z — a;) c«/ 2 a + (x — y) cy z {3 + (y — z) cy 2 y) 
— (cy*a + c*/ 2 /3+ cy % y)((z — x) cya -\- (x — y) cyfi + (y — z) cyy) , 

which reduces to (2y — x — z) a + (2z — x — y) (3 -\- (2x — y — z) y . Hence 

i S*.v'p. V (v', cy''p — p) + S.tfp. V(cy>*p, p) = — -L (3tf.p 2 — S 2 .v'p) tf.v'p. 

Dividing out the common factor v', equation (36) reduces to 

£ ftp = — i [ F(p, «') cos Q-g- + 4>^) + F(«/, cy 2 p) cos <?>] + -g- v'S.v'p, 

*Cf. Tait's Quaternions, Art. 66. 



or 
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3(up = v'S.v'p — 2 [ F(p, v') cos (~y + <?>) + V(v', cy'*p) cos <£>] , (41) 



a formula which is very similar to (22), to which it reduces when v' = v. 

When <£> = 90° and p is perpendicular to i/, equation (38) becomes 

3 W =V3F(p, v'), (42) 

as it should. 

Taking S.v'p = 0, we have now 

3cy'*y = V(vf, p) cos (-j- + &) + V{cy' % 9 , v>) cos 3; (43) 

and differentiating this with respect to S-, we find that 

^.3cy^p = — 7(t/, p) sin (-^ 4- &) d$ — F(cy 2 p, t/) sin » 

3 ( s +4) . (44) 

27t 

8* 

Writing for brevity cy n " — <f, so that qp 2w = cy' s p = p, and observing that 
F(«', p) = xg~2p, and V(cy'*p, v')=V(q s p, #') = — q~*q3p=z — g«p we can write 

equations (39) and (40) in the form 

_ — / 2th \ 5? r 

Sq*p =q 2 p COS ( -y- + S) — § • p cos 3 (45) 

and <*.<fp = / +f -P^- (46) 

We now perceive that taking a a unit vector perpendicular to v', the equation of 

a plane curve may be written 

p = rq^a (47) 

where r is a scalar. Differentiating twice in succession with respect to the time, 
we find that 

p' = riq»a + rq* + h.^-, (48 ) 

p"=(''-<t)V+( 2 '! +^§y +f - <*•> 

These equations give the velocity and acceleration of a body moving in a plane 
resolved along and perpendicular to the radius vector.* 

* Cf. Tait's Quat., Art. 337. 
48 



378 Chapman : On the Matrix which Represents a Vector. 

Having thus, by a succession of opei^ations comprised well within the domain 
of algebra, and with the aid of notations borrowed from Quaternions, obtained 
an operator which rotates a vector through a definite angle about a definite axis, 

I will now pause to compare q & with the unit versor e w of quaternions.* 

?* 
e 17 and q & have precisely the same effect on a vector in the plane perpen- 

dicular to their axes ; but s" operating on g gives e * , a versor whose angle 

exceeds that of e™ by — - , its plane being the same ; while q & operating on its 

axis v' simply leaves it unchanged. For this reason the effects of these opera- 
tors on vectors not perpendicular to their axes are different, that of q & being the 
simpler. We may also compare q & with the quaternion operator q.a,q~ x which 
rotates ot about the axis and through double the angle of q.f 
Returning to equation (41) and observing that if 

S = aa, + b{3 + cry, h' = da. + V{3 + c'y, (50) 

then V ( S, S>) = (be' — cb') a + {pa! — ad) {3 + (aV — ba') y , 

so that 

V(p + 2cy'*p, v') = (2x — y — z)a + (2y — z — x) + (2a — y — x)cc; 

and writing the equation 

3^p = v'S.v'p + cos 4> V (p + 2cy'*p , v') + V3 sin <p V (p , v) , (51) 

we have, if $ is a very small angle, and if we note that v'.Sv'p = (x -f- y -f- z) 

(a + + y), 

3^p= 3p + V3^7(p, v'), (52) 

to quantities of the second order. 

For a rotation of <r through an angle dS about an axis v", which meets v r , 
we have 3f j a _ 3<y + ^3^7(5, v"). 

In this writing jup for a , we shall have the effect of the two rotations in succes- 
sion upon p. This gives 

V/wp = 3 (3p + ^/3dq>V(p, v') + +/3d&V(3p + */3dcpV(p, v 1 ), v") j 

= 9p + 3*/3dc}>V(p, «/) + 3^3c^F(p, OJ j" ( 53 ) 

= 9^*^ d *p. ' 

Here the order of the rotations is indifferent, but that is of course not true for 
finite rotations, ft' being in general a function of (ip. 

*Taifc, p. 88 ; Ex. 13. t Tait, p. 261. t Compare Tait, p. 260. 
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4. — Applications to Trigonometry. 



The determinant 



cya cya cy*a 



cy(3 cy(3 cy*(3 
cyy cyy cy\ 

is identically zero. Hence we conclude that for any three vectors whatever 
oyaVtf, y) + cy@V(y, a) + cyyV(a, (3) = 0. 

If a, (3, y lie in a plane and their lengths are a, b, c and the angles between them 
A, BondC, then F(/?, y) = bcv*siiaA, V(y, a) = cav' sin /?, V(a, /?) = a6w / sin G, 
«' being the direction of the perpendicular on the plane, and equation (50) 
becomes, after dividing by «/, 

be sin A cya + ca sin Bcyfi + ab sin Ocyy =■ 0; 

operating on this with cy % , we have, after dividing by abc, and denoting the ver- 
sors* of a, (3, y by a', (3', y', 

a' sin A + [3' sin B + y' sin (7=0, (55) 

which shows that if the sides are proportional to the sines of the opposite angles, 

the plane figure will be closed and a triangle. Ifa + {3 + y=0 correspond to 

any other closed figure formed by multiples of a', (3', y', it must be a consequence 

of (51); otherwise we could eliminate either vector and must conclude that the 

remaining two are paralleif 

Again, from the equation 

-a = P + y, (56) 

we have, by squaring and taking scalars, 

S.a 2 =S.p*+ Sy+ 2S.(3y, (57) 

or a? = b % +c 2 —2bccosA. (58) 

The use of the negative sign in the last term is necessary from the manner of 
measuring the angle A. 

With a spherical triangle we may proceed as follows: We have, if a, (3, y 
are unit vectors to the vertices, 

S. V (a, (3) V(J3, y) = — cos B sin a sin c. 



*Tait, p. 33. t Cf. Gibbs, Vector Analysis, p. 49. 
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Now 

V(a, /3) V (/3, y) = (cyacy 2 (3 — cy{3cy*a)(cy{3cy 2 y — cyycy 2 (S) 

= — cyf&cy 2 (ay) -f- cy (@y) cy 2 (a/3) + cy (a/3) c# 8 (/3y) — cy (a/) c?/ 2 /? 3 ; 

and recalling that 

S.cy (a/3) cy 2 (/3y) = S.fiycy* (a/?) and tf.c?/ (/3^) cy 2 (a/3) = S.(3ycy (a/3), 

by a formula of Art. 2, we may write this 

S.V(a, /3) F(/3, y ) = S.py(cy(aP) + c^ 8 (a/3)) — ^./3 2 (cy(a r ) + ^(a^)) 

= /S^y (vS.a(3 — a/3) — #./? 2 (v#.ay — ay) 
= #. ^7^. a/3 — S. (PS. ay 
= cos a cos c — cos b, 

by previous formulae. We conclude that 

cos b = cos a cos c + sin a sin c cos B.* 

This extremely simple process may be compared with the analysis given in Tait, 
p. 56, Art. 90, which must be read to get the result on p. 71. 

Without following so closely the familiar processes of quaternions, we may 
proceed as follows : 



V(a, /3)F(/3, y)=z 



cya cy(3 ^ cy(3 cyy 
cy 2 a cy*(3 cy 2 (3 cy 2 y 

cy (a/3) + cy 2 (a/3) , cy (ay) + cy 2 (ay) 
cy(3 2 + cy 2 ? 2 , cy (0y) + cy 2 (py) 

vS.afi — a/3, vS.ay — ay 

vS.{3* — {3 2 , vS.(3y—(3y 

= vS.apS.(3y — (3yS.a(3 — afiS.py — vS.^S.ay 
+ ayS.p? + (3 2 jS.ay. 

Hence, taking scalars, we have 

S.a(3S.(3y — S.^S.ay — S.V(a, /3) F(/3, y) 

as before, by remembering that S.v = 3. 



*Tait, p. 71. 



